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Abstract. We study the distribution of baryonic and lu- 
minous matter within the framework of a hierarchical sce- 
nario. Using an analytical model for structure formation 
which has already been checked against observations for 
galaxies, Lyman-a clouds, clusters and reionization pro- 
cesses, we present its predictions for the bias of these ob- 
jects. We describe its dependence on the luminosity (for 
galaxies or quasars) or the column density (for Lyman-a 
absorbers) of the considered objects. We also study its red- 
shift evolution, which can exhibit an intricate behaviour. 
These astrophysical objects do not trace the dark matter 
density field, the Lyman-a forest clouds being undercor- 
related and the bright galaxies overcorrelated, while the 
intermediate class of Lyman-limit systems is seen to sam- 
ple the matter field quite well. We also present the dis- 
tribution of baryonic matter over these various objects. 
We show that light does not trace baryonic mass, since 
bright galaxies which contain most of the stars only form 
a small fraction of the mass associated with virialized and 
cooled halos. We consider two cosmologies: a critical den- 
sity universe and an open universe. In both cases, our 
results agree with observations and show that hierarchi- 
cal scenarios provide a good model for structure forma- 
tion and can describe a wide range of objects which spans 
at least the seven orders of magnitude in mass for which 
data exist. More detailed observations, in particular of the 
clustering evolution of galaxies, will constrain the astro- 
physical models involved. 



Key words: cosmology: large-scale structure of Universe - 
galaxies: evolution - quasars: general - quasars: absorption 
lines - intergalactic medium - galaxies: clusters 



1. Introduction 

In the standard cosmological scenario large-scale struc- 
tures observed in the present universe were formed by 
the amplification through gravitational instability of small 
primordial density fluctuations (usually assumed to be 
gaussian). Within hierarchical models (as in the CDM 



case: Peebles 1982) small scales collapse first to form 
bound objects which later merge to build more massive 
halos as larger scales become non-linear. This describes 
the evolution under the action of gravity of the dark mat- 
ter component. However, observed objects (e.g., quasars, 
galaxies, Lyman-a clouds...) are formed by baryonic mat- 
ter (gas, stars) which does not necessarily follow the be- 
haviour of the dark matter component since it undergoes 
additional processes (e.g., radiative coohng). Moreover, by 
looking at specific astrophysical objects (e.g. bright galax- 
ies) one selects peculiar dark matter environments where 
the former are most likely to arise. Thus, to relate obser- 
vations to the underlying dark matter component (which 
is usually gravitationally dominant) and to test scenarios 
for structure formation one needs to study the link be- 
tween the dark and baryonic components. This can also 
constrain the astrophysical models used to build galaxies 
and other objects from the density field. 

One tool to describe the distribution of objects in the 
universe, in addition to their abundance, is their two-point 
correlation function which measures their clustering. Then 
one often defines the bias of these halos as the ratio 
of this quantity to the dark matter correlation function, 
which shows in a clear quantitative way whether these 
objects are more or less clustered than the dark matter 
component. In this article, we present the bias (and the 
correlation length) associated with various objects (galax- 
ies, quasars, Lyman-a clouds, clusters). The formulation 
of this bias is obtained (Schaeffer 1985,1987; Bernardeau 
& Schaeffer 1992,1999) within the framework of a general 
description of structure formation in the universe provided 
by a scaling model (Schaeffer 1984; Balian & Schaeffer 
1989) which describes the density field in the non-linear 
regime. These predictions specifically take into account 
the deep non-linearity of the density field and are concep- 
tually different from those obtained using linear theory 
(e.g.. Mo & White 1996). Indeed, the latter authors eval- 
uate the bias of primordial overdensities that will collapse 
later, using linear theory, with the implicit assumption 
that this bias is conserved during the course of evolution. 
They find the bias depends separately on the mass, ra- 
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dius and overdensity. Here, we use the bias calculated di- 
rectly using non-linear theory. This bias is seen to depend 
on very specific internal properties of the considered non- 
linear objects, namely on a unique parameter x, quite close 
to the internal velocity dispersion within the objects, thus 
exhibiting a scaling law. Such a specific dependence has 
indeed been found to hold in numerical simulations (Mun- 
shi et al. 1999b), and we consider this approach to provide 
a more trustworthy description of the present clustering 
properties as well as their evolution at higher redshifts. 

In the present paper, we specifically study the de- 
pendence of the bias on the observable properties of the 
considered objects (luminosity, column density) and its 
evolution with redshift. In contrast with most previous 
works, our analytical model relating dark to observable 
mass is not an ad-hoc parameterization. Indeed, it is sim- 
ply a consequence of a general model (Valageas & Scha- 
effer 1997) which has already been applied to galaxies 
(Valageas & Schaeffer 1999), Lyman-a clouds (Valageas 
et al. 1999), clusters (Valageas & Schaeffer 2000) and 
reionization studies (Valageas & Silk 1999a, b). Thus, its 
predictions have already been compared with many ob- 
servations (e.g., galaxy luminosity function, column den- 
sity distribution of Lyman-a clouds, amplitude of the UV 
background radiation field, cluster X-ray luminosity func- 
tion,...). Hence there is no free parameter chosen in this 
article. Next, in addition to the bias we also present the 
distribution of baryonic matter over the various objects 
we describe. 

This article is organized as follows. In Sect|| we de- 
scribe our prescription for mass functions, while in Sect.^ 
we recall how the bias is obtained in this model and we 
compare our formalism with other prescriptions. In Sect.|| 
we present our results for the case of a critical density uni- 
verse. We first describe the clustering of galaxies, quasars, 
Lyman-a clouds and clusters and then we study the repar- 
tition of baryonic matter within various classes of objects. 
Finally, in Sect.|^ we present the case of an open universe. 



Jeans length. In any case, we write the multiplicity func- 
tion of these objects as (see Valageas & Schaeffer 1997 for 
details): 



M M X 



(1) 



Here, the parameter x associated with a halo of mass M 
at redshift z is defined by: 



1 + A(M, z) 
x(Ai, z) = = 

^ ' m{M,z),z] 



(2) 



where 



' 6(ri,r2) with V=-7tR^ 



is the average of the two-body correlation function 
C2(ri,r2) over a spherical cell of radius R and provides 
the measure of the density fiuctuations in such a cell (thus 
typical objects have a; ~ 1). The scaling function H{x) de- 
pends on the initial power-spectrum of the density fiuctua- 
tions and must be taken from numerical simulations. How- 
ever, from theoretical considerations (Balian & Schaeffer 
1989) it is expected to satisfy the asymptotic behaviour: 



a; < 1 : H{x) cx x'^ 

with Lo ~ 0.5, LOs ^ 
it must satisfy 



x H{x) dx — 1 



a; > 1 : H{x) oc x'" 



-3/2, ^ 10 to 20, and by definition 



(3) 



These properties have been checked for various 
P(fc) by Bouchet et al (1991), Colombi et al. 
(1992,1994,1995,1997), Munshi et al. (1999a), while 
the mass functions obtained from (|^) for various con- 
straints A(M) (for a constant A which was taken from 
—0.5 to 5000 or for (1 -I- A) (x M) have been shown to 
provide reasonable approximations to numerical results 
in Valageas et al. (2000). 



2. Multiplicity functions 

Since we wish to describe the properties of various classes 
of objects like Lyman-a clouds and galaxies which are de- 
fined by specific constraints we need a formalism which 
can handle more general multiplicity functions than the 
usual mass function of "just-virialized halos" . To this or- 
der, we shall assume that the non-linear density field obeys 
the scaling model detailed in Balian & Schaeffer (1989). 
Thus, we define each class of objects (Lyman-a clouds or 
galaxies) by a constraint A(M, z) on the density contrast 
of the underlying dark-matter halo of mass M at redshift 
z. The "just-virialized halos" correspond to the special 
case A(M, z) = lS.c{z) with Ac ~ 177. For our purposes 
we consider the cases of a constant density contrast, aris- 
ing from a virialization constraint, and of a constant radius 
(i.e. (1 -|- A) cx M), arising from a fixed cooling radius or 



3. Bias 

3.1. Our formulation 

Within the framework of the scale- invariance of the many- 
body matter correlation functions ^p(ri, r^) which led 
to the mass function (||) , one can show that in the highly 
non-linear regime the bias characteristic of two objects 
factorizes and is a function of the sole parameter x intro- 
duced in Sect.|[ see Bernardeau & Schaeffer (1992,1999). 
Thus we write the correlation function of objects of type 
(1) with objects of type (2) as: 

ei,2(r) = &(xi) b{x2) (4) 
with 

X < 1 : 6(x) (X x^^"'^^/^ and x > 1 : 6(x) cx x . (5) 
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This behaviour has been sucessfuUy checked in numerical 
simulations by Munshi et al. (1999c). This allows us to 
obtain the redshift evolution of the bias associated with 
galaxies, Lyman-a clouds or quasars, once we define the 
constraints A(M, z) which characterize these various as- 
trophysical objects (that in turn define the proper value 
of X to be used). 

3.2. Comparison with some other models 

We can note that alternative models have been proposed 
in the litterature to describe the clustering of galaxies. 
Here we briefly compare our prescription to such mod- 
els, in order to clarify their main differences. Thus, semi- 
analytic models of galaxy and dark matter clustering 
based on Press-Schechter-like ideas have recently been 
proposed (Seljak 2000; Peacock & Smith 2000). They de- 
scribe the dark matter density field as a collection of 
smooth halos with a mean density profile. The mass func- 
tion of these "just-virialized" halos (defined by an over- 
all density contrast A — Ac ^ 177) is obtained from a 
slightly modified Press-Schechter mass function (Press & 
Schechter 1974). By construction, such a model neglects 
the substructures of virialized halos since they are de- 
scribed by an average density profile. This is not a problem 
at large scales where the correlations among galaxies are 
governed by the distribution of their host virialized halos 
(note that since the Press-Schechter prescription is based 
on the linear density field these correlations should actu- 
ally apply to very large scales which are still in the linear 
regime). However, at scales of the order of the size of large 
"just-virialized" halos the clustering of galaxies is set by 
the distribution of these galaxies within those larger ob- 
jects (which actually correspond to clusters of galaxies at 
large masses). 

Hence one needs to add a specific prescription to cope 
with this regime. Thus, one introduces the mean number 
{N){M) of galaxies contained in halos of mass M. This 
function is derived from N-body simulations or observa- 
tions. Then, one puts a galaxy at the center of the halo 
and the {N — 1) remaining galaxies at random within the 
halo, following the mean dark matter density profile. In 
order to include the fluctuations of N{M) one can even 
reproduce the second moment {N{N — 1))(M) obtained 
from simulations (Seljak 2000). However, it is clear that 
the predictive power of such a model is quite limited. In- 
deed, the output of the model (i.e. the two-point correla- 
tion function of galaxies) is almost directly related to the 
input (the first two moments of the galaxy distribution 
within virialized halos). Thus, this can at best serve as a 
test of N-body simulations if this data is taken from sim- 
ulations (while it would not make much sense to obtain 
this data from the observations one would later on try to 
recover). In particular, in order to derive the p— point cor- 
relation function of galaxies in this line one would need 
the moment (A^^) of order p... On the other hand, if one 



only includes the mean {N){M) (taken from observations 
for instance) as input, one can argue that the model tests 
whether the spatial galaxy distribution is well described by 
the mean dark matter density profile. However, as shown 
in Valageas (1999) such a model contradicts the results of 
numerical simulations for the dark matter distribution it- 
self when one takes into account the constraints provided 
by the first five order moments of counts-in-cells statis- 
tics. This implies that one could not get higher-order mo- 
ments (e.g., the three-point galaxy correlation function) 
from this model. 

Moreover, it is not clear how much information one 
can draw from the mere comparison of the galaxy power- 
spectrum predicted by the model with observations. In- 
deed, any sensible model should provide a reasonable 
agreement with the data since it should more or less follow 
the dark matter power-spectrum which is itself normal- 
ized to observations of large-scale clustering (through ag). 
Then, it seems difficult to obtain narrow constraints on 
the input parameters from only one curve, in view of the 
number of these free parameters and the approximations 
involved in such models (e.g., importance of the scatter of 
galaxy properties, shape of the distribution P{N) within 
virialized halos, role of dark matter substructures,...). 

Such approaches differ from ours in many respects. 
First, we use a mass function derived from a scaling 
ansatz for the dark matter density field (Balian & Scha- 
effer 1989) which applies to the non-linear regime. It is 
fully specified by the sole function H{x), obtained from 
counts-in-cells statistics, and by the behaviour of the non- 
linear power-spectrum. Most importantly, it allows us to 
get various multiplicity functions for different classes of 
objects in addition to "just-virialized" halos, defined by 
a relation A(M) or R{M) which specifies their overall 
density contrast or radius as a function of their mass 
(Valageas & Schaeffer 1997; Valageas et al. 2000). This 
enables us to bypass the link between galaxies and "just- 
virialized" halos which hinders the approaches described 
above. Indeed, we directly recognize the galaxies them- 
selves, defined by a relation R{M) obtained from cooling 
constraints (Valageas & Schaeffer 1999, 2000). 

Then, as explained above the clustering of these ob- 
jects is obtained (Bernardeau & Schaeffer 1992,1999) from 
the behaviour of the p— point correlation functions 
which led to the mass function (|l|). Hence it does not re- 
quire any additional information. Moreover, although we 
restrict ourselves to the two-point correlation functions, 
our formulation also predicts the p— point correlation func- 
tions of any order p for all the objects we study here 
(galaxies, QSOs,...). This clearly illustrates the predictive 
power of our model. Besides, although we do not take into 
account the scatter of galaxy properties with respect to 
their mass, our procedure automatically includes the part 
of the scatter of the galaxy distribution with respect to 
the distribution of "just-virialized" halos which is due to 
the statistics of the dark matter itself. 
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Finally, the spirit of our present study is rather differ- 
ent from the works which are mainly intended to describe 
galaxy clustering (Seljak 2000; Peacock & Smith 2000) as 
we merely compare with observations the predictions of 
a general model of structure formation which was fully 
constrained in previous papers which dealt with its conse- 
quences for galaxies, QSOs, Lyman-a clouds and clusters. 
In particular, we have no free parameters left to improve 
the agreement with the observed galaxy clustering and we 
check our model against a vast array of observations, rang- 
ing from Lyman-a clouds up to clusters, through galaxies 
and QSOs. Thus, our goal is not to use observations of 
galaxy clustering to constrain some free parameters, but 
to check whether a sensible scenario of structure forma- 
tion can match this data. Besides, since we use an analytic 
model where the various physical processes can be clearly 
identified, we can get some valuable insight into the clus- 
tering properties of the range of objects we study here. For 
instance, we shall see that a model like ours can explain 
in simple physical terms the qualitative and quantitative 
differences between the redshift evolution of the cluster- 
ing of bright galaxies and faint galaxies, small Lyman— a 
forest clouds and damped systems, or QSOs and galaxies. 

Another method to study the clustering of galaxies is 
to use N-body simulations, which may be supplemented 
by semi-analytic modelling of star- formation or galaxy for- 
mation (e.g., Benson et al. 2000a). In principle, this allows 
one to avoid the assumptions and approximations involved 
in an analytic model like ours to describe the dark mat- 
ter density field. However, this technique may be limited 
by numerical resolution and computer time constraints. In 
any case, this is certainly a promising method but it does 
not provide yet the clustering properties of the full range 
of objects we study here in a unified fashion. 

4. Critical density universe 

We first consider the case of a critical density universe 
rJm = 1 with a CDM power-spectrum (Davis et al. 1985), 
normalized to as — 0.5. We choose a baryonic density 
parameter ftf, — 0.04 and Ho ~ 60 km/s/Mpc. We use 
the scaling function obtained by Bouchet et al. (1991) for 
H{x). 

4-.1. Galaxies 

4.1.1. Present universe 

We use the model of galaxy formation and evolution de- 
scribed in detail in Valageas & Schaeffer (1999), where it 
was checked against many observations (luminosity func- 
tion, TuUy-Fisher relation...). We define galaxies by two 
constraints: 1) a virialization condition A > Ac where 
Ac ^ 177 is given by the usual spherical model, and 2) 
a cooling constraint tcooi < which requires the gas to 
cool within a few Hubble times at formation. At high virial 



temperature and low redshift, the condition 2) is the most 
stringent and can be approximated by a constant radius 
constraint R = Rcooi- At high redshift it becomes irrel- 
evant since any halo which satisfies 1) also satisfies 2). 
These two conditions define for galaxies the constraint 
A(M, z) we introduced in Sect.|[ Once we have identified 
in this manner the dark-matter halos which correspond 
to galaxies we still need to set up a model for star for- 
mation to obtain the stellar properties of these objects. 
As detailed in Valageas & Schaeffer (1999) our model in- 
volves 4 components: (1) short-lived stars which are re- 
cycled, (2) long-lived stars which are not recycled, (3) a 
central gaseous component which is deplenished by star 
formation and ejection by supernovae winds, replenished 
by infall from (4), a diffuse gaseous component. This al- 
lows us to obtain the luminosity of these galaxies. Since 
observations usually refer to the mean bias (or correlation 
function) of galaxies above a given luminosity, we define: 

^ /r..(i)f ' 

where we used the factorization of the bias, see (||). We 
present in Fig.|l| the dependence of the bias on the B-band 
luminosity. 
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Fig. 1. The luminosity dependence of the bias 6(> Lb) at 
z = 0. The solid line is the mean bias b{> LB)/b{> L*) 
of galaxies brighter than Lb normalized to the bias of 
galaxies brighter than Af* = —19.5. The data points are 
from Benoist et al. (1996) (crosses) and Willmer et al. 
(1998) (squares). 

We see that we recover the trend of larger bias for 
brighter galaxies, although our dependence is somewhat 
steeper than the observed relation. However, we note that 
in our analytical model each luminosity is associated with 
a specific object of a given mass, radius,... In a more real- 
istic model one should include some scatter in this relation 
which would lessen the final luminosity dependence of the 
bias (for a large scatter such that the luminosity is uncor- 
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related with the properties of the underlying dark matter 
halo we would have no B-magnitude dependence left). 

From the relation (^) we can also obtain the cross- 
correlation of galaxies of luminosity Lb with all galaxies, 
in a fashion similar to . We can convert the bias b evalu- 
ated in this manner into a correlation length ro{b) defined 
by: 

ro(6) = 6^/^ ro , (7) 

where 7 = 1.8 is the approximate slope of the correlation 
function and tq is the correlation length of the dark matter 
density field. Indeed in our model the correlation functions 
of various objects and of the matter distribution all have 
the same slope in the non-linear regime (note that Dave 
et al. 1999 find from numerical simulations that the bias 
shows very little scale dependence, if any, in the non-linear 
regime). This prescription allows us to compare in Fig.^ 
our results to observations performed by Loveday et al. 
(1995). 
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Fig. 2. From left to right the correlation length tq at z ~ 
of the correlation function of galaxies between the magni- 
tude bins -22.8 < Mb < -20.8, -20.8 < Mb < -19.8 
and —19.8 < Mb < —15.8 with all galaxies (i.e. with 
-22.8 < Mb < -15.8). The point on the far right is 
the correlation length of all galaxies with —22.8 < Mb < 

— 15.8. The data points (crosses) are from Loveday et al. 
(1995) and the squares are our predictions. 

We see that our results are compatible with the obser- 
vations. We find an increase of the bias (hence of tq) for 
brighter galaxies in accordance with the data, except for 
the last point at —22.8 < Mb < —20.8 where the bias is 
measured to be slightly smaller than for —20.8 < AIb < 

— 19.8. Nevertheless, since this contradicts the trend ob- 
served by Benoist et al. (1996) this disagreement is prob- 
ably not conclusive. 

We present in FigJ^ the correlation function of galax- 
ies brighter than Mb = —18.6 at z = 0. Our results agree 
reasonably well with observations, although the slope we 
obtain may be a bit too steep at large scales. This is a 
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Fig. 3. The correlation function of galaxies brighter than 
Mb — —18.6 at z = (solid line). The dashed line is the 
linear prediction. The data points are from Baugh (1996) 
(squares) for —21.6 < Mb < —18.6 and from Loveday et 
al. (1995) (crosses) for -22.8 < Mb < -15.8. 

well-known discrepancy for CDM initial conditions with 
rim = 1 and could be improved with a slightly different 
initial power-spectrum. The dashed line is the correlation 
function given by the linear theory multiplied by the rele- 
vant bias factor. Hence it coincides with the actual correla- 
tion function at large scales. At small scales the non-linear 
dynamics increases the clustering as compared to the lin- 
ear predictions, in agreement with observation. Note that 
in principle the non-linear predictions are valid for ^ ^ 1, 
but are in practice seen to hold as soon as ^ > 1. 

4.1.2. Redshift evolution 

We now turn to our predictions for the redshift evolution 
of the bias associated with galaxies recognized by their 
B-band luminosity, shown in Fig.H. 




-22 -20 -18 -16 



Fig. 4. The redshift evolution of the relation luminosity- 
bias. The solid line corresponds to z = 0, the dashed line 
to z = 2 and the dot-dashed line to z = 4. 
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At "low" luminosities Mb > —21, the bias grows with 
redshift since galactic halos correspond to increasingly 
rare overdense regions (because structure formation was 
not as developed as it is today). For very bright galaxies 
Mb < —21 we recover the same trend except at low red- 
shifts where the bias suddenly increases with time. This 
evolution is due to a luminosity effect. Indeed, in our 
model at low z < 1 very massive galaxies see their lu- 
minosity decrease with time as their initial gas content 
gets exhausted by star formation (which is very efficient 
in this high-density environment). As a consequence, when 
selecting galaxies by a fixed B-band luminosity, one looks 
at more massive, higher virial temperature halos at z = 
than at z = 1. This shift to higher density fluctuations 
acts opposite to the usual trend seen for lower luminosity 
galaxies and leads to the increase with time of the bias at 
low z for these very bright magnitudes. 
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Fig. 5. The redshift evolution of the bias associated with 
galaxies below various B-band magnitude thresholds. The 
solid line corresponds to Mb < —20, the dashed line to 
AIb < —18 and the dot-dashed line to Mb < —15. 

We present directly in Fig.|^ the redshift evolution of 
the bias associated with galaxies brighter than three B- 
band magnitude thresholds. This clearly shows the de- 
pendence on luminosity of the redshift evolution itself. 

We show in Fig|| the redshift evolution of the comov- 
ing correlation length rg associated with galaxies brighter 
than various magnitude thresholds. The two data points 
at z < 0.1 correspond to z = but they have been sHghtly 
shifted to be more clearly seen. We can see that our re- 
sults agree reasonably well with observations, except for 
the data points from Le Fevre et al. (1996) which how- 
ever are inconsistent with data from Carlberg et al. (1998) 
and Postman et al. (1998). The comoving number den- 
sity of Lyman-break galaxies at z = 3 is n ^ 2.7 x 10"'^ 
Mpc~^ (see Adelberger et al. 1998) which is also the 
abundance of galaxies brighter than Mb ~ —18.7 (and 
M > 5 X 10" Mq) in our model. The square at z = 3 
shows our prediction for the clustering strength of galax- 
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Fig. 6. The redshift evolution of the correlation length tq 
(in comoving Mpc) associated with galaxies below various 
B-band magnitude thresholds. The solid line corresponds 
to Mb < -20, the dashed line to Mb < -18 and the dot- 
dashed line to Mb < —15. The data points are from Le 
Fevre et al. (1996) (crosses, -21.3 < Mb < -18.3), Carl- 
berg et al. (1998) (filled squares, Mb < —19.5), Postman 
et al. (1998) (circles, / < 20) and Maghocchetti & Mad- 
dox (1999) (open triangles). The filled triangle at z = 3 
corresponds to Lyman-break galaxies (observations from 
Adelberger et al. 1998) while the open square is our pre- 
diction for galaxies which have the same abundance as 
Lyman-break galaxies. 

ies with this abundance. We see that this simple identi- 
fication leads to good agreement with observations (filled 
triangle) . The comoving correlation length of galaxies first 
decreases at higher z, following the decline of the dark- 
matter two-point correlation function. This trend is even 
larger for very bright galaxies since it is amplified by the 
decline of their bias, as shown in Fig.^. However, at high 
redshifts z > 2, increases with z because of the strong 
growth of the bias. Of course, this trend is stronger for 
brighter galaxies. We note that contrary to Kauffmann et 
al. (1999) we obtain a "dip" in the evolution of tq, while 
these authors only found this feature for a low-density flat 
cosmology. This is explained by the slower growth with z 
of the bias associated with galaxies in our model, as com- 
pared to theirs, which is insufhcient to override at low z 
the decrease of clustering measured by the decline of the 
dark-matter correlation function. Thus, such a feature de- 
pends on the astrophysical model used to build galaxies. 
However, we note that observations show that decreases 
with larger z in the range < z < 0.5, in agreement with 
our model. We recall here that our results are direct conse- 
quences of a detailed model for structure formation which 
has already been compared with many observations, so 
that no attempt was made to "improve" the agreement 
of our predictions with the data. In fact, as shown by the 
discrepancy between various observations, the latter still 
have a significant uncertainty. However, we recover the 
general trend of a decrease of tq up to z ^ 2 and a sub- 
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sequent growth at higher redshifts. This behaviour also 
agrees with the results obtained by Baugh et al. (1999) 
from a semi- analytic model. 

4-2. Quasars 

From the description of galaxies presented above we have 
also developed a model for the quasar luminosity function, 
as detailed in Valageas & Silk (1999a). We assume that a 
fraction Ag < 1 of galactic halos host a quasar with a mass 
Mq proportional to the mass of gas Mgc available in the 
inner parts of the galaxy: Mq = F Mgc. For most galaxies 
this also implies Mq ~ F Mg where Mg is the stellar mass 
and we use F = 0.01. We also assume that quasars shine 
at the Eddington limit with a radiative efficiency e — 0.1, 
which means that the quasar life-time is tQ = 4.4 e 10* 
yr. This determines the quasar bolomctric luminosity: 



Ln = 



eMc 



(8) 



We present in Fig.^ the redshift evolution of the quasar 
luminosity - bias relation. Here b{> Lb,z) is the bias of 
quasars brighter than Lb sX redshift z. 



o O 



z = 



log(Lg/L^) 



Fig. 7. The redshift evolution of the bias associated with 
quasars above various luminosity thresholds. The solid line 
corresponds to z = 0, the dashed line to z = 2 and the 
dot-dashed line to z = 4. 

We see that the bias increases with redshift since at 
higher z quasars of a given luminosity correspond to rarer 
density fluctuations. The high luminosity cutoff seen in the 
figure comes from the fact that in our model very massive 
galaxies have already consumed most of their gas content 
at low z, as explained in the previous section, so that the 
maximum quasar luminosity declines with time because 
of fuel exhaustion. Of course, at higher z the quasar lu- 
minosity function cutoff also declines because of the knee 
of the mass function of collapsed objects (Valageas & Silk 
1999a). 

From the bias b{> Lb,z) we can derive the comov- 
ing correlation length of quasars. This is shown in Fig.H. 



1.5 



0.5 




Fig. 8. The redshift evolution of the comoving correlation 
length associated with quasars above various luminosity 
thresholds. The dashed line corresponds to Lb > 10^^ Lq 
and the solid hue to Lb > 10^^ Lq. The two crosses are 
observations from La Franca et al. (1998) for 0.3 < z < 
1.4 (lower point) and 1.4 < z < 2.2 (upper point). The 
triangle shows the result at 2.7 < z < 4.75 of observations 
by Stephens et al. (1997). 

We can check that our results agree reasonably well with 
observations. In particular, we recover the increase with 
redshift of the correlation length. La Franca et al. (1998) 
noticed that this behaviour differs from the decrease of rp 
observed for galaxies from z = to z = 1 which we also ob- 
tain, as shown in Fig.^. Thus a simple hierarchical model 
like ours is able to simultaneously explain both of these 
redshift dependences. The stronger increase with redshift 
of the correlation length associated with quasars, as com- 
pared with galaxies, also translates into the faster growth 
of their bias factor, see FigJJ and Fig.|^. This difference 
between galaxies and quasars is due to the redshift evolu- 
tion of their respective ratios of (mass)/ (luminosity). As 
described above, in our model the luminosity of a quasar 
is proportional to its mass of gas, see (^), and to the mass 
of stars in the host galaxy. Hence the bias of quasars of a 
given luminosity grows strongly with redshift since objects 
with a fixed mass of central gas are increasingly rare at 
earlier times when structure formation is less advanced. 
This overrides the decrease of the dark matter correla- 
tion length with z and it leads to a rise of ro(> Lb), as 
seen in Fig.^. On the other hand, galaxies with a con- 
stant mass of star-forming gas had a higher luminosity in 
the past because the ratio Mg/L decreases with redshift 
(see for instance Fig. 17 in Valageas & Schaeffer 1999). 
This is due to the fact that the smaller age of the uni- 
verse tn selects more massive and brighter stars with a 
life-time oc tn. In other words, at later times the con- 
tribution from faint long-lived stars which progressively 
accumulated in the galaxy ever since it was born becomes 
increasingly important as compared with the contribution 
of bright newly-born short-lived stars which quickly dis- 
appear as supernovae and thus cannot accumulate. This 
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implies that by looking at a fixed luminosity one selects a 
smaller mass of gas and stars at higher redshift. 

Thus, if we start at z = with galaxy and quasar 
luminosities which correspond to the same host halos, at 
large redshift z > the galaxies with the same luminosity 
correspond to smaller halos than the quasar hosts. Hence 
galaxies correspond to less extreme objects and have a 
smaller bias. Thus at low z their correlation length follows 
the decreases of the dark matter tq, until z 1, beyond 
which the rise of their bias becomes sufficiently strong to 
lead to a growth of ro(> Lg), see Fig.^. The precise ampli- 
tude of these redshift evolutions depends somewhat on the 
astrophysical models used for quasars and stars (IMF,...). 
However, the agreement we obtain with these observations 
is remarkable since our model was built independently of 
these considerations, with no free parameter left for this 
present discussion. Thus this behaviour appears as a nat- 
ural outcome of such models which suggests they provide 
a realistic description of the global properties of star and 
quasar formation, even though the details of these pro- 
cesses are still poorly known. 

4-3. Lyman— a absorbers 

As shown in Valageas et al. (1999), the formalism intro- 
duced in Sect. II also allows us to build a model for Lyman- 
a clouds, some of which are objects with a lower density 
than the mean density of the universe. We describe these 
absorbers as three different populations of objects. Low 
density halos with a small virial temperature see their 
baryonic density fiuctuations erased over scales of the or- 
der of the Jeans length Rj. These mass condensations 
form our first class of objects, defined by the scale Rj, 
which can be identified with the Lyman-a forest. They 
correspond to low density clouds which can even be un- 
derdense (but stick out within very low-density regions or 
voids). On the other hand, potential wells with a large 
virial temperature maintain their density profile. Thus, 
they constitute a second population of absorbers where 
the column density observed along a line of sight depends 
strongly on the impact parameter. They can be identified 
with the Lyman-limit systems. Finally, the deep neutral 
cores of these halos (because of self-shielding) form our 
third class of objects which correspond to the damped 
systems. 

We present in Fig.^ the redshift evolution of the bias 
associated with Lyman-a absorbers above various column 
density thresholds. 

We note that large column density absorbers have a 
bias greater than unity since they correspond to high den- 
sity fluctuations (massive halos) while forest clouds which 
are low density fluctuations located in underdense areas 
have a small bias 6 < 1. For Nhi < lO^"* cm~^ the bias de- 
creases at larger redshift because the higher mean density 
of the universe implies that to observe a constant small 
column density one must look at increasingly underdense 
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16 
log(N„,) 



Fig. 9. The redshift evolution of the bias associated with 
Lyman-a absorbers above various column density thresh- 
olds. The solid line corresponds to z = 0, the dashed line 
to z = 2.5 and the dot-dashed line to z = 5. 

objects. In a similar fashion, for Lyman-limit or damped 
systems the bias first declines from z = Otoz — 2.5. How- 
ever, since these objects correspond to virialized halos the 
bias eventually grows with z as the structure formation 
process is less advanced. Of course there is also an inter- 
mediate range of Nhi where the behaviour is more intri- 
cate. We display the redshift evolution of the bias asso- 
ciated with clouds above three different column densities 
thresholds in FigJlCj. 



1.5 



N„, 10'-" c:rii- = 



N„, > 1 .6 1 O" orri- = 



2 3 



Fig. 10. The redshift evolution of the bias associated with 
Lyman-a absorbers above various column density thresh- 
olds. The solid line corresponds to Nhi > lO^"* cm~^, the 
dashed line to Nhi > 1-6 10^^ cm~^ and the dot-dashed 
line to Nhi > 2 lO^o cm-^. 



We see that the redshift evolution of the bias of Lyman- 
a absorbers can display a non-trivial behaviour, since 
many processes play a role in the properties of these 
clouds, including the amplitude of the UV background 
radiation field. Thus, a phenomenological parameterisa- 
tion of the form b (x N'^j (1 -I- z)^ is unlikely to provide 
an accurate description. Note however that intermediate 
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clouds in the range 10 < Nhi < 10 cm have a bias 4.5. Distribution of matter 



which remains reasonably close to unity at all redshifts 
z < 5 and thus can be used as tracers of the matter den- 
sity field. A comparison of our results with observations 
for the two-point correlation function in velocity space is 
described in Valageas et al. (1999). 

4-. 4. Clusters 



T a 5 kcV 



2 — 



T > 1 keV 



T a 0.5 keV 



The models we have built in previous studies for galax- 
ies (Valageas & Schaeffer 1999) and Lyman-a clouds 
(Valageas et al. 1999) allow us to obtain the redshift evo- 
lution of the fraction of baryonic matter associated with 
various components. 




Fig. 11. The redshift evolution of the comoving correla- 
tion length associated with clusters above various virial 
temperature thresholds. The dot-dash line corresponds to 
T > 0.5 keV, the dashed line to T > 1 keV and the solid 
line to T > 5 keV. The two squares are observations from 
Croft et al. (1997) at 2; < 0.18 for clusters with a richness 
7^ > 50 (lower point) and 7?. > 110 (upper point). The 
triangle shows the result at z < 0.2 of observations by 
Borgani et al. (1999) for the X-ray brightest Abell cluster 
sample. 

Finally, we show in Fig.|l^ the redshift evolution of 
the correlation length associated with clusters above vari- 
ous virial temperature thresholds. Here we simply defined 
clusters as "just-virialized" halos characterized by the 
usual density contrast Ac ~ 177. A detailed description 
of our model for clusters and the evolution of their X-ray 
luminosity function is presented in Valageas & Schaeffer 
(2000). We see that our results agree reasonably well with 
observations. As was the case for quasars, the rise with 
redshift of the bias associated to these halos leads to an 
increase of their correlation length. However, this increase 
is slower than for quasars, see Fig.^, because here we select 
these halos by a constant temperature threshold. This im- 
plies that their mass evolves as M (x (1 -I- z)~^/^ hence we 
consider smaller halos at high z which are more common 
(and have a smaller bias) than objects of constant mass. 
Thus, we have obtained a simple model of structure forma- 
tion which is consistent with observations from Lyman-a 
clouds (M ~ 10*^ Mq) up to clusters (M ~ 10^^ Mq), i.e. 
over 7 orders of magnitude. 



Fig. 12. The redshift evolution of the fraction of bary- 
onic mass embedded within Lyman-a forest clouds (up- 
per shaded area), galactic halos (white area) and stars 
(lower shaded area). The dotted lines separate the contri- 
butions from bright {Mb < —19.5, lower part) and faint 
(Mb > —19.5, upper part) galaxies for stars and gas, re- 
spectively. 



We present in Fig.|T^the contributions of stars, galactic 
gas and Lyman-a forest clouds. In our model the sum of 
these three constituents is unity. Indeed, all the matter 
is embedded within mass condensations which we identify 
with galaxies (deep potential wells where the gas can cool) 
or with Lyman-a forest clouds (small density fiuctuations 
heated by the UV background radiation) . We see that the 
mass fraction within stars is always small ('^ 10% at z = 
0) while the contribution of Lyman-a forest clouds shows 
a fast increase with redshift to constitute ~ 70% of the 
matter at z = 5 when there are few deep potential wells. 
Bright galaxies {Mb < —19.5) contain most of the mass 
embedded within stars but they only constitute a small 
part of the gas which was able to cool. Moreover, their 
total mass of gas starts declining at low redshifts z < 1 
as it is converted into stars. This trend does not appear 
for faint galaxies because their star formation process is 
less efficient while they accrete gas which was previously 
part of the Lyman-a forest clouds. Moreover, some of these 
small galaxies have just formed as new potential wells have 
built up. Note that if the IGM is reheated by supernovae 
or quasars up to T ~ 5 10^ K, which is suggested by 
the observed relation temperature - X-ray luminosity of 
clusters (e.g., Ponman et al. 1999), the fraction of matter 
which has been able to cool and appears as galactic gas 
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Table 1. Baryonic parameter fli associated with various 
objects at z = and z = 3, for a critical universe with 
fib — 0.04. The second line which appears for some con- 
tributions is the observational estimate from Fukugita et 
al. (1998). 





z = 


2 = 3 


stars 


0.0057 


0.0006 




0.002 - 0.0057 




galactic gas 


0.024 


0.015 


groups 


0.012 


0.00002 




0.011 - 0.045 




clusters 


0.0034 


0.00002 




0.0014 - 0.0044 




Lyman- a forest 


0.01 


0.025 






0.01 - 0.05 



in Fig.|2| would be smaller. This was studied in details in 
Valageas & Silk (1999b). 




13 3 4 



Fig. 13. The redshift evolution of the fraction F of bary- 
onic mass embedded within galaxies (upper solid line: to- 
tal mass, lower solid line: gas in the core), groups (dashed 
line), clusters (dot-dashed line) and stars (dotted line). 

We show in Figjl^ the mass fraction we associate with 
groups and clusters of galaxies, as well as with stars and 
galaxies. We define groups as virialized halos (i.e. with 
a density contrast Ac) which do not satisfy the con- 
straint tcooi < tn- They correspond to massive potential 
wells with a large virial temperature (at high T we have 
tcooi ~ VT)- At high redshifts this contribution declines 
since an increasing proportion of virialized halos can cool 
(indeed tcooi o<: 1/p while tn oc l/i/p), see also Valageas 
& Schaeffer (1999). We identify clusters as groups with a 
virial temperature T > 1 keV. At low z they only form a 
small sub-class of galaxy groups while at high z all groups 
have large virial temperatures T > 1 keV as explained 
above. Note that in our model the mass associated with 
groups is also associated to the member galaxies, so that it 
is a part of the fraction of matter recognized as embedded 



within galaxies, shown in Fig. ^ and in FigjT^ (upper solid 
line). A detailed study of groups and clusters in presented 
in Valageas & Schaeffer (2000). 

From the fraction of matter Fi associated with vari- 
ous objects and the baryonic density parameter Qh — 0.04 
we used in our model, we also obtain the baryonic den- 
sity parameters fti = FiUf, for various classes of objects. 
We display in Tabj^ a comparison of our results with ob- 
servational estimates by Fukugita et al. (1998). We note 
that our predictions are consistent with the data. In Tabj^ 
we removed from the mass embedded within clusters and 
groups shown in Figjl^ the fraction of baryons in stars or 
galactic cores which they contain. 

5. Open universe 

Using the same formalism we can also study the case of 
an open universe fim — 0.3, J^a = 0, with a CDM power 
spectrum normalized to erg = 0.77. We use fif, = 0.03 and 
Hn = 60 km/s/Mpc. 

5.1. Galaxies 




— 1 o 1 

log(R) [ Mpc ] 



Fig. 14. The correlation function of galaxies brighter than 
Mb — —18.6 at z = as in FigJ^. The upper solid line 
corresponds to the open universe with Q,^ = 0.3 and f2A = 
0, while the lower dashed-line is the linear prediction. The 
upper dashed line shows the low-density flat universe with 
rim = 0.3 and = 0.7. The middle solid line (resp. 
dashed hne) displays the case Om = 0.5 and = (resp. 
= 0.5 and Qa = 0.5). 

We show in Fig.|l^ the correlation function of galax- 
ies brighter than AIb — —18.6 at z = (upper solid 
line). Our predictions are similar to the results obtained 
in Figj^. However, a bump appears at intermediate scales 
[R ^ 2 Mpc) where ^ ~ 10. This corresponds to the sharp 
increase of £,{R) from the linear regime (^ ^ 1) to the 
non- linear regime (^ ^ 100). This leads to a "step- like" 
profile which is also seen in numerical simulations (e.g.. 
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Valageas et al. 2000) and which gets larger for low den- 
sity universes (Peacock & Dodds 1996). This is due to the 
slow-down of the linear growth factor while the two-point 
correlation function keeps increasing as ^(i?,a) oc in 
the highly non- linear regime (Valageas & Schaeffer 1997). 
The observed two-point correlation function on the other 
hand is surprisingly close to a power-law as soon as ^ > 1. 
Such a behaviour requires a precise balance between the 
curvature of the initial power-spectrum, the distortions in- 
duced by the non-linear dynamics and the bias of galaxies 
(which also depends on non-gravitational processes which 
affect galaxy formation). Note however that using numer- 
ical simulations Pearce et al. (1999) and Benson et al. 
(2000a) obtain a good match to observations as the galax- 
ies they get are antibiased relative to the dark matter 
on small scales which "conspires" to provide a power-law 
galaxy correlation function. This behaviour is partly due 
to the spatial exclusion of halos. Moreover, the expression 
(^ we use for the bias is only valid at large separations 
r ^ i? between objects of size R. Thus, the discrepancy at 
small scales might be cured by a more rigorous calculation 
which would not use the approximation r 3> i?. 

On the other hand, note that in the case of a criti- 
cal density universe this "bump" in the mildly non-linear 
regime is much smaller, which leads to good agreement 
with observations (see Fig.^). In particular, we also dis- 
play in Fig.^ the dark matter correlation function we 
would obtain for a low-density flat universe (upper dashed 
line) with 0^ = 0.3 and SIa = 0.7, with a slightly dif- 
ferent power-spectrum from Bardeen et al. (1989) (and 
r = 0.18) using the non-linear fit from Peacock & Dodds 
(1996). Thus this corresponds to a bias of unity, which is 
a good approximation for galaxies with Mb < —18.6 at 
z = 0, see FigJ^ and Fig. 15, The main point is that the 
"bump" is much weaker for a flat cosmology, as can also be 
seen in Peacock & Dodds (1996). The intermediate lines 
show the dark matter correlation functions obtained for 
^^m = 0.5 (with f2A = and ^Ia = 0.5) and normalized 
to cTg ~ 0.6 (and F = 0.3). Thus, the shape itself of the 
galaxy correlation function does not allow one to strongly 
discriminate between these various cosmological scenarios. 
Although the "bump" which appears for the open uni- 
verse with rijii = 0.3 seems to disfavour this model, the 
inaccuracy involved in the modelling of the galaxy clus- 
tering properties is certainly too large to draw meaningful 
conclusions, especially since numerical simulations suggest 
that simple physical processes could remove this feature 
(Benson et al. 2000a). 

We present in Fig.^ the redshift evolution of the 
luminosity-bias relation. We can see that our results are 
again similar to those obtained in FigJ^, since the astro- 
physical model is the same and it must satisfy the same 
observational constraints at low z. However, the bias of 
bright galaxies increases more slowly with redshift than 
for a critical density universe (compare the curve at z = 4 
with Fig.^. This is due to the higher normalization as 
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Fig. 15. The redshift evolution of the relation luminosity- 
bias. The solid line corresponds to z = 0, the dashed line 
to z = 2 and the dot-dashed line to z = 4. 

of the power-spectrum and to the slower redshift evolu- 
tion of the linear growth factor which means that massive 
galactic halos correspond to less extreme density fluctua- 
tions at z = 4 than for the case ri„i = 1 when they have a 
similar clustering pattern at z = 0. 
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Fig. 16. The redshift evolution of the correlation length rp 
(in comoving Mpc) associated with galaxies below various 
B-band magnitude thresholds, as in Fig.^ for an open 
universe. 



We show in Fig.^the redshift evolution of the comov- 
ing correlation length ro associated with galaxies below 
various B-band magnitude thresholds. The galaxy cor- 
relation length is larger than for the previous case of a 
critical density universe because the normalization of the 
power-spectrum (erg) is slightly higher. This also implies a 
slower redshift evolution. We note that the clustering pat- 
tern we predict for Lyman break galaxies is somewhat too 
strong as compared to observations. However, this could 
be due to some extinction by dust of the luminosity of 
these LBG. Indeed, in our model the number of galaxies 
brighter than Mb = —18 is n = 3 x 10~^ Mpc~^ while 
the abundance of LBG is n ~ 8.5 x lO"'^ Mpc'^ (Adel- 
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berger at al. 1998). This means that we would recover the 
right correlation length if we assumed that LBG corre- 
spond to 30% of these galaxies with Mb < —18. An alter- 
native would be that LBG correspond to galaxies which 
are in the midst of a starburst phase (e.g., Somerville et 
al. 2000) which increases their luminosity (note that the 
absolute LBG magnitude is rather large: Mn < —21.8 for 
an apparent magnitude TZ < 25.5 without K-correction). 
In this case, Lyman-break galaxies would also be drawn 
from a larger parent population which could again lead to 
a smaller correlation length. Note on the other hand that 
it would be difficult to reconcile the model with observa- 
tions if we had predicted a smaller value of tq than the 
one given by the data. 




Fig. 17. The redshift evolution of the fraction of bary- 
onic mass embedded within Lyman-a forest clouds (up- 
per shaded area), galactic halos (white area) and stars 
(lower shaded area). The dotted lines separate the contri- 
butions from bright {Mb < —19.5, lower part) and faint 
{Mb > —19.5, upper part) galaxies for stars and gas. 
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Fig. 18. The redshift evolution of the fraction F of bary- 
onic mass embedded within galaxies (upper solid line: to- 
tal mass, lower solid line: cores), groups (dashed line), 
clusters (dot-dashed line) and stars (dotted line). 

Table 2. Baryonic parameter Qi associated with various 
objects at 2 — and z — 3, for an open universe with 
= 0.3, Qb — 0.03. The second line which appears 
for some contributions is the observational estimate from 
Fukugita et al. (1998). 





2 = 


2 = 3 


stars 


0.0039 


0.0006 




0.002 - 0.0057 




galactic gas 


0.019 


0.013 


groups 


0.011 


0.0007 




0.011 - 0.045 




clusters 


0.003 


0.00002 




0.0012 - 0.0044 




Lyman- a forest 


0.007 


0.016 






0.01 - 0.05 



6. Conclusion 



5.2. Distribution of matter 



We show in Fig. 17 the contributions of stars, galactic gas 
and Lyman-a clouds to the baryonic density in the uni- 
verse as a function of redshift. We obtain a behaviour sim- 
ilar to Fig.|l|. 

We display in Fig.|8| the redshift evolution of the frac- 
tion of matter embedded within galaxies, groups, clusters 
and stars. We see that the mass fraction formed by groups 
and clusters declines more slowly than in Fig.^. Again 
this is due to the higher normalization trg and to the slower 
redshift evolution of the linear growth factor. 

We show in Tab.^ the baryonic density parameters Hi 
associated with various objects. We note that our results 
still agree reasonably well with observations. 



In this article we have described the redshift evolution 
of the clustering of various objects (galaxies, quasars, 
Lyman-a clouds, clusters) within the framework of a uni- 
fied model of structure formation which has already been 
checked against many observations. Hence there are no 
free parameters chosen to fit the observations considered 
in this article, since the model was already overconstrained 
by previous studies. We have shown that the bias of these 
astrophysical objects displays an intricate behaviour which 
cannot be described by simple power-law parameteriza- 
tions. Indeed, the link between these objects and the un- 
derlying dark matter component reflects the influence of 
many factors (cooling processes, star formation efficiency, 
evolution of the UV background radiation field...) which 
can have competing effects. We have shown that in our 
model the bias associated with galaxies strongly increases 
with their luminosity. This was already emphasized in the 
earlier calculations (Schaeffer 1987; Bernardeau & Scha- 
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effer 1992) as being a genuine outcome of the hierarchical 
clustering hypothesis. We recover the trend seen in ob- 
servations, although our dependence may be slightly too 
strong. However, this could be due to the simplicity of our 
model which does not include any scatter in the properties 
of the objects we model. The redshift evolution we obtain 
also agrees with data, in particular the correlation length 
of Lyman-break galaxies observed at z = 3 is consistent 
with our results. 

We have also described the bias associated with 
Lyman-a clouds. Thus, it appears that absorbers in the 
range 10^^ < Nhi < 10^^ cm~^ have a bias close to unity 
(hence they trace the dark matter density field) in contrast 
to other clouds and galaxies which show a very different 
behaviour of their clustering properties. Our results also 
agree with observations for quasars and clusters. In par- 
ticular, they explain in a natural fashion the qualitative 
and quantitative differences between the redshift evolu- 
tion of the clustering of QSOs and of galaxies. Thus, our 
model provides a good description of the intrinsic proper- 
ties and the spatial clustering of a wide variety of objects, 
from small underdense Lyman-a absorbers up to massive 
clusters, which covers six orders of magnitude in mass. 

Finally, we have described the redshift evolution of the 
distribution of matter over the various objects observed 
in the universe. Our results show a reasonable agreement 
with observations. This is not surprising since our model 
has already been checked against the galaxy luminosity 
function, the number density of Lyman-a clouds, which 
constrain the amount of matter embedded within these 
objects. In particular, it appears that light does not trace 
baryonic mass as bright galaxies {Mb < —19.5) which 
contain most of the stars only provide a small fraction of 
the mass associated with virialized and cooled halos. 

We have also performed the same analysis for the case 
of an open universe. In both cosmologies our results agree 
reasonably well with observations but the details of as- 
trophysical models are still too uncertain to discriminate 
between these two cases. However, for the open universe 
the galaxy correlation function we obtain shows a shoulder 
which does not appear in observations. This could be due 
to the approximation involved in our calculation which is 
only valid on scales much larger than the size of the ob- 
jects. In particular, some numerical simulations do not get 
this "bump" (Benson et al. 2000a). On the other hand, 
they do not recover the luminosity dependence of the 
galaxy bias (Benson et al. 2000b) and simulations cannot 
simultaneously describe the clustering of the whole range 
of objects we studied in this paper. Nevertheless, the fact 
that our predictions are nearly consistent with the data 
over such a wide range of objects shows that hierarchi- 
cal scenarios provide a very good description of structure 
formation. On the other hand, more detailed observations 
of the amplitude of the clustering of galaxies selected by 
their luminosity will constrain the astrophysical models. 
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